Oi  a 

00  I 

—  w 

Oi  g 

rH 


Bethatda,  Maryland  20084 


DAVID  W.  TAYLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER 


DESCRIPTION  OF  COMPUTER  PROGRAM  CAB3DYN  FOR 
THE  TIME  DOMAIN  ANALYSIS  OF  THE  THREE-DIMENSIONAL 
LONG-TIME  DEPLOYMENT  BEHAVIOR  OF  GENERAL  OCEAN 
CABLE  SYSTEMS 


Henry  T,  Wang 


Raymond  S.  Cheng 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


SHIP  PERFORMANCE  DEPARTMENT 


DEPARTMENTAL  REPORT 


C  p£  H 
H  a:  V) 
H  H 
X  if. 
M  UJ 
X  X  W 
H  X 
it*  cO 

X  X 
X  G  V 


December  1981 


DTNSRDC/SPD-06 33-03 


NOW  OTNSRDC  5602  30  2  80) 
i supersedes  3960  46) 


8  2 


MAJOR  DTNSRDC  ORGANIZATIONAL  COMPONENTS 


NOW  DTNSRDC  3960/43  IR«*  2  f 

‘  ^  4.  ***M**W«a*VS*Sg»n^  "  J 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  PAGE  <WI ton  0«t»  E nlcr.d) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2,  GOVT  ACCESSION  NO.I  !  RECIPIENT'S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


report  number 


SPD-0633-03 


4  TITLE  (and  Submit) 

DESCRIPTION  OF  COMPUTER  PROGRAM  CAB3DYN  FOR  THE 
TIME  DOMAIN  ANALYSIS  OF  THE  THREE-DIMENSIONAL  LONG¬ 
TIME  DEPLOYMENT  BEHAVIOR  OF  GENERAL  OCEAN  CABLE  I  performing  oro.  report  number 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 


FINAL 


7.  AUTHORfiJ 


I.  CONTRACT  OR  GRANT  NUMBCRfaj 


Henry  T.  Wang  and  Raymond  S0  Cheng 


1  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

David  W.  Taylor  Naval  Ship  Research  and 
Development  Center 
Bethesda,  Maryland  20084 
_  C0NTR0LUNG  orf[Ct  name  and  address  ~~  " 

Naval  Air  Development  Center  (Code  3043) 
Warminster,  Pennsylvania  18974 


12.  REPORT  DATE 

December  1981 


u.  number  of  pages 


MONITORING  AGENCY  NAME  A  AOORESS  (It  dllloront  from  Controlling  Olllct)  IS.  SECURITY  CLASS,  (ol  thlo  report) 

UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


[T«  DISTRIBUTION  STATEMENT  (ol  IM*  Rapori) 


Approved  for  Public  Release:  Distribution  Unlimited 


17  DISTRIBUTION  STATEMENT  (of  tho  mbttroct  ontorod  In  Block  20,  ti  dlfioront  from  Report) 


I  II  SUPPLEMENTARY  NOTES 


II-  KEY  WORDS  (Contlnum  on  roveteo  «Jd«  If  neceeeery  md  Identity  by  block  number) 

Buoy-Cable-Body 
Thr  ee-D iuen  s iona 1 
Time  Domain 
Deployment  Behavior 

20  ABSTRACT  (Continue  on  reveree  el  do  II  neceeeaty  and  Identity  by  block  nuaiSarJ 

*nThe  report  describes  in  detail  Program  CAB3DYN,  which  analyzes  in  the 
time  domain  the  three-dimensional  long-time  deployment  behavior  of  free- 
floating  and  towing  ocean  cable  systems.  The  report  describes  the  coordinate 
systems  used,  the  modeling  of  the  time-varying  current  profile,  and  the 
formulations  for  the  surface  buoy  or  prescribed  surface  motion,  intermediate 
bodies,  and  the  connecting  cable.  Detailed  input  instructions  include  a 
listing  of  the  input  READ  statements,  the  definition  of  each  input  variable  , _ 


00  ^ 


EDITION  OF  1  NOV  AS  IS  OBSOLETE 

S-'N  0102-LF-OU-6601 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dal*  InlaradJ 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whan  Data  EntararfJ 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OFTHIS  PAGEfWian  Data  Enlatad; 


TABLE  OF  CONTENTS 


PAGE 


ABSTRACT  -  1 

ADMINISTRATIVE  INFORMATION  -  1 

INTRODUCTION  -  I 

CALCULATION  OF  CABLE  TENSION  -  4 

PREVIOUS  APPROACHES  -  4 

PRESENT  APPROACHES  -  6 

COORDINATE  SYSTEMS  -  7 

CURRENT  PROFILE - - -  3 

PRESCRIBED  SURFACE  MOTIONS  -  9 

SURFACE  BUOY  FORCES -  10 

INTERMEDIATE  BODY  FORCES - • -  12 

CABLE  FORCES -  14 

FINAL  EQUATIONS  OF  MOTION -  19 

DESCRIPTION  OF  COMPUTER  PROGRAM  -  20 

MAIN  PROGRAM -  20 

SUBROUTINE  DYNA -  21 

SUBROUTINE  CUR -  21 

SUBROUTINE  KUTMER  -  22 

SUBROUTINE  XLINEAR  -  22 

SUBROUTINE  PRESC  -  22 

PROGRAM  STORAGE  AND  TIME  REQUIREMENTS -  22 

INPUT  INSTRUCTIONS -  23 

READ  STATEMENTS -  23 

DEFINITION  OF  INPUT  VARIABLES  -  24 


iii 


PAGE 


DEFINITION  OF  INPUT  VARIABLES  FOR  SURFACE  BUOY, 

(FSM  (1)  >  2000) -  27 

EXPLANATORY  NOTES  -  27 

SAMPLE  PROBLEMS  - 28 

FREE-FLOATING  BUOY-CABLE  PROBLEM  -  28 

TOWING  CABLE  PROBLEM  -  34 

SUMMARY -  36 

ACKNOWLEDGMENT  -  38 

REFERENCES -  39 

LIST  OF  FIGURES 

1  -  Geometrical  Configuration  of  Cable  System - - - — — - -  41 

2  -  Segmented  Representation  of  Cable - * -  42 

3  -  Parameters  for  Free-Floating  Buoy-Cable  System  of 

Sample  Problem  1  - - - - - — - - - - -  43 

4  -  Parameters  for  Cable  Towing  Lower  Array  of  Sample 

Problem  2 -  44 

LIST  OF  TABLES 

1  -  Final  Steady-State  Results  for  Free-Floating 

Buoy-Cable  Syetem  -  45 

« 

2  -  Values  of  Velocity  x  in  Ft/Sec  Calculated  by  Program 

CABMOD  During  Deployment  of  Free-Floating  Buoy-Cable 

System - 46 

3  -  Values  of  Cable  Angle  <{>v  in  Deg  During  Deployment  of 

Cable  Towing  Lower  Array  at  1  Knot  - - * -  47 

iv 


PAGE 


Steady-Stata  Values  of  Cable  Angle  <j>  In  Deg 

n  i 1  .  i J  i n _ _ _ /iininniu  ^ i r  _  J:  x  


Calculated  by  Program  CABTRAN  for  Various 

Values  of  V  and  NCAB - - - - - - -  48 

S 

Values  of  Cable  Angle  (f>  in  Deg  During  Deployment 

of  Cable  Towing  No  Arra^  at  1  Knot - - - - - - - - -  50 

Values  of  Cable  Angle  <p  in  Deg  During  Deployment 

of  Cable  Towing  No  Area}1  at  5  Knots - - — - - -  52 


v 


ABSTRACT 


The  report  describes  In  detail  Program  CAB3DYN,  which  analyzes 
in  the  time  domain  the  three-dimensional  long-time  deployment  behavior 
of  free-floating  and  towed  ocean  cable  systems.  The  report  describes  the 
coordinate  systems  used,  the  modeling  of  the  time-varying  current  profile, 
and  the  formulations  for  the  surface  buoy  or  prescribed  surface  motion, 
intermediate  bodies,  and  the  connecting  cable.  Detailed  input  instructions 
include  a  listing  of  the  input  READ  statements,  the  definition  of  each  input 
variable,  and  comments  on  usage.  Two  sample  problems  on  the  deployment  of 
free-floating  and  towed  cable  systems  are  given  to  illustrate  usage  of  the 
program,  calculated  results,  and  computer  execution  time. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  jointly  supported  by  the  Naval  Air 
Development  Center,  Code  3043,  under  Work  Request  N62269/WR/00726  dated, J17. __ , 
July  1980  and  by  the  Naval  Sea  Systems  Command,  Code  03R,  under  the  Ship  Per¬ 
formance  and  Hydromechanics  Program.  The  work  was  performed  under  internal 
work  units  1-1552-548,  1-1552-533,  and  1-1507-101. 

INTRODUCTION 

This  report  presents  a  detailed  description  of  Program  CAB3DYN,  which 
analyzes  in  the  time  domain  the  three-dimensional  long-time  deployment  behavior 
of  free-floating  and  towed  ocean  cable  systems.  The  term  long-time  dynamic 
behavior  refers  to  motions  where  time  scales  are  large  enough  so  that  inertia 
effects  are  not  of  major  importance.  Examples  of  such  motions  are  the 
deployment  of  a  free-floating  cable  system  from  an  arbitrary  initial  state  to 
the  final  steady-state  and  the  behavior  of  a  towing  cable  undergoing  a  prescribed 
maneuver. 


1 


To  some  extent,  the  formulation  in  Program  CAB3DYN  may  be  considered  to  be 

an  extension  of  the  two-dimensional  dynamic  formulation  contained  in  Program 
1  2 

CABMOD  *  to  three  dimensions  by  using  the  coordinate  system  of  the  three- 

3  4 

dimensional  steady-state  Program  CAB3E  *  .  The  cable  system  consists  of  a 
buoy  or  towing  platform  at  the  surface,  a  cable  whose  properties  may  vary  along 
its  length,  and  the  presence  of  intermediate  bodies  along  the  cable .  The  current 
profile  is  modeled  by  a  series  of  input  magnitudea  and  directions  at  various 
depths,  as  in  Program  CAB3E.  The  cable  is  represented  by  a  series  of  straight 
segments,  as  in  Program  CABMOD. 

However,  there  are  significant  differences.  Those  factors  which  are  pertur¬ 
bations  to  the  overall  mean  long-time  motion  are  not  considered  In  the  present 
program.  In  particular ,  ocean  waves  and  rotational  modes  of  the  surface  buoy 
are  neglected.  These  result  in  considerable  simplification  of  the  formulation 
and  reduction  in  the  number  of  input  variables  for  the  surface  buoy.  On  the  other 
hand,  hydrostatic  buoyancy  and  drag  characteristics  of  the  buoy  are  modeled  in 
somewhat  greater  detail.  Another  considerable  simplification  results  from  the 
neglect  of  added  inertia  for  the  cable.  This  eliminates  coupling  between  the 
equations  of  motion  for  a  given  cable  node.  For  two  sample  problems  it  is  shown 
that  the  cable  added  inertia  terms  have  nearly  negigible  effect  on  the  motion. 

The  report  starts  by  giving  a  brief  review  of  cable  studies  which  are  most 
closely  related  to  the  present  study.  The  report  then  describes  the  coordinate 
systems  used  in  the  program,  the  modeling  of  the  time-varying  current  profile, 
and  the  formulations  for  the  surface  buoy  (or  prescribed  surface  motion),  inter¬ 
mediate  bodies,  and  the  connecting  cable.  A  description  is  given  of  the  various 
subroutines  of  the  program.  Detailed  input  instructions  are  given,  including  a 
listing  of  the  input  READ  statements,  the  definition  of  each  input  variable,  and 
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comment 8  on  the  entry  of  input  data.  Wherever  possible,  the  input  variables  are 
matched  to  those  contained  in  Program  CABMOD. 

Two  sample  problems  are  given  to  illustrate  usage  of  the  program,  output, 
and  computer  execution  time.  The  two  problems  involve  the  deployment  of  a  free- 
floating  buoy-cable  system  and  a  towing  cable  from  an  initially  vertical  position 
to  the  steady-state  configuration  in  the  presence  of  various  current  profiles 
and  towing  speeds,  respectively.  For  each  problem,  results  computed  by  Program 
CAB3DYN  are  compared  to  those  obtained  by  Program  CABMOD  for  two-dimensional  cases. 
It  is  shown  that  when  the  option  of  neglecting  cable  added  inertia  terms  In  Program 
CABMOD  is  exercised,  the  two  programs  agree  Identically  for  the  towing  cases  and 
agree  in  an  average  sense  for  the  free-floating  cases.  The  presence  of  the  buoy 
pitch  terms  in  Program  CABMOD  introduces  a  perturbation  to  the  overall  motion 
which  is  not  modeled  by  Program  CAB3DYN.  For  both  problems,  it  is  further  shown 
that  there  Is  little  difference  between  the  motions  predicted  by  Program  CABMOD 
with  and  without  cable  added  Inertia  effects. 

For  the  free-floating  buoy-cable  problem,  an  efficient  approach  is  indicated 
for  using  Program  CAB3DYN,  in  conjunction  with  Program  CAB3E,  to  predict  the 
detailed  configuration  of  the  system  in  the  presence  of  a  slowly  varying  current 
profile.  For  the  towing  cable  problem,  results  calculated  by  Programs  CAB3DYN 
and  CABMOD  are  compared  to  those  calculated  by  Program  CABTRAN  ,  which  is  a  two- 
dimensional  program  which  completely  neglects  Inertia  effects.  It  is  shown  that 
for  cases  where  the  cable  is  towing  a  neutrally  buoyant  lower  array,  Programs 
CABMOD  and  CAB3DYN  predict  unstable  motions  of  the  lower  segment  at  higher  towing 
speeds  while  Program  CABTRAN  predicts  incorrect  steady-state  configurations.  For 
cases  of  a  cable  alone,  with  no  lower  array,  the  results  predicted  by  Programs 


3 


CABMOD  and  CAB3DYN  agree  reasonably  well  with  those  predicted  by  Program  CABTRAN. 

The  report  concludes  by  giving  a  brief  summary  of  the  results  for  the  two 
sample  problems. 

CALCULATION  OF  CABLE  TENSION 

PREVIOUS  APPROACHES 

It  is  well  known  that  for  most  cable  systems,  computer  time  and  coat 

are  determined  by  the  need  to  track  the  high-frequency  longitudinal  vibrational 

mode  along  the  cable.  In  particular,  the  maximum  time  step  Atu  which  may 

M 

be  taken  is  given  by 

AiM  Jp. l/c'  to 

where  p  is  the  mass  per  unit  length  of  the  cable 

L  is  the  length  of  the  cable  segment 
Cj (■  AE)  is  the  elastic  constant  of  the  cable 

A  is  the  cross-sectional  area  of  the  cable 

E  is  the  elastic  modulus 

For  transient  ot  other  high-frequency  motlunB,  it  is  of  interest  to  carefully 
model  the  longitudinal  vibrations.  However,  for  the  present  case  of  deployment 
motions,  whose  time  scales  are  orders  of  magnitude  longer  than  that  of  the  longi¬ 
tudinal  vibrations,  it  is  of  interest  to  avoid  the  time-step  restriction  imposed  by 

Equation  (i).  Thus,  previous  studies  which  consider  this  aspect  were  carefully 
investigated. 
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For  the  most  part,  these  studies  attempt  to  increase  the  allowable  time  step 

by  using  the  condition  of  inextensibility  to  obtain  the  tension,  or  avoid  the 

tension  calculation  altogether  by  considering  only  motions  normal  to  the  cable. 

6 

Walton  and  Polachek  add  to  the  equations  of  motion  for  the  cable  additional 

nonlinear  equations  on  the  displacement  which  express  the  condition  of 

inextensibility.  At  each  time  step,  the  tensions  along  the  cabls  are  iterated 

a  number  of  times  until  the  inextensibility  condition  is  satisfied.  Froide«aux 
7 

and  Scholten  also  use  an  iterative  procedure  to  calculate  the  tension.  Their 
procedure  consists  of  continually  guessing  the  tension  at  the  top  of  the  cable 
and  calculating  the  tension  along  the  cable  by  taking  the  longitudinal  motion 
to  be  the  same  along  the  cable  until  the  calculated  tension  below  the  last  body 
is  nearly  zero. 

g 

Ketchman  and  Lou  show  that  for  cases  where  the  surface  motion  is  prescribed 

and  the  lower  body  is  cable-like,  such  as  an  array,  the  tension  can  be  directly 

calculated  along  the  entire  cable,  starting  from  the  lower  body.  However,  for 

cases  where  the  lower  body,  or  any  body  along  the  cable,  is  not  cable-like  (i.e., 

has  drag  and  added  inertia  properties  which  depend  on  fixed  spatial  directions 

rather  than  on  cable  inclination)  the  calculation  of  tension  would  again  not  be 

straightforward.  Similar  complexities  would  arise  if  the  surface  motion  is  not 

prescribed,  as  in  the  caBe  of  a  surface  buoy. 

9  5 

Rupe  and  Thresher  and  Paul  and  Soler  ,  use  methods  which  solve  for  only  the 
motions  normal  to  the  cable  and  thus  avoid  the  need  to  calculate  the  tension.  After 
the  normal  motions  have  been  calculated,  the  tension  may  be  obtained,  if  desired. 
Rupe  and  Thresher  consider  the  cable  to  be  a  compound  pendulum  composed  of  a  series 
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of  links.  They  solve  for  the  angular  coordinates  by  using  Lagrange's  Method. 

Paul  and  Soler  use  a  unique  approach.  They  show  that  for  typical  towing  maneuvers , 
it  is  permissible  to  neglect  all  Inertia  effects,  thus  reducing  the  equations  of 
motion  to  first  order.  They  then  show  that  it  is  possible  to  start  at  the  bottom 
link  and  use  a  quasi-static  approach  which  balances  drag  and  weight  forces  to 
obtain  the  normal  velocities  along  the  entire  cable.  However,  as  in  the  case  of 
Reference  8,  both  of  the  above  approaches  work  well  for  only  the  case  where  the 
surface  motion  is  prescribed  and  the  cable  is  free  of  bodies  or  has  only  cable¬ 
like  bodies.  The  presence  of  non-cable-like  bodies  would  make  their  formulations 
more  complex  and  the  presence  f  a  surface  buoy,  whose  motion  is  not  prescribed, 
would  cake  the  calculation  of  tension  an  essential  part  of  the  procedure. 


PRESENT  APPROACH 

In  view  of  the  Iteration  procedure  required  in  References  6  and  7,  and  the 
lack  of  generality  of  the  approaches  given  in  References  5,  8,  and  9,  it  was 
decided  to  take  the  tension  T  to  be  proportional  to  the  strain  £  ,  as  follows: 

T=  Tr  +  ( 


where  T  is  the  reference  tension  at  which  the  strain  is  taken  to  be  zero, 
r 

In  Program  CABMOD,  there  is  an  Internal  damping  term  proportional  to  strain 
rate  e  ,  and  e  is  raised  to  an  arbitrary  power  C^.  Computer  cost  may  be 
minimized  by  using  the  reduced  modulus  technique*®  Introduced  by  the  first 
author.  The  technique  essentially  maximizes  the  allowable  time  step  At^ 
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given  In  Equation  (1)  by  replacing  the  actual  elastic  constant  by  the  lowest 
permissible  value  which  will  still  give  accurate  results.  For  deployment  pro¬ 
blems,  it  appears  reasonable  to  use  values  of  C  which  will  lead  to  maximum 
strains  of  0.01.  A  one  percent  elongation  should  have  little  effect  on  the 
transverse  deployment  motion. 

COORDINATE  SYSTEMS 


The  two  coordinate  systems  used  in  the  computer  program  are  shown  in  Figure 
1.  As  mentioned  in  the  Introduction,  these  coordinate  systems  are  identical  to 

4 

those  used  In  Program  CAB3E  . 

The  equations  of  motion  for  the  cable  system  are  derived  for  the  spatial 
(x,  y,  z)  coordinate  system  whose  origin  lies  at  the  undisturbed  ocean  surface. 

The  x  and  z  axes  lie  in  the  horizontal  plane  and  the  y  axis  is  positive  In  the 
direction  of  gravity.  Cable  drag  Is  most  conveniently  derived  for  a  second 
coordinate  system  attached  to  the  cable,  (X,  Y,  Z) ,  where  the  Y-axis  is  directed 
along  the  cable.  The  cable  coordinate  system  is  obtained  by  first  rotating  the 
spatial  system  by  the  angle  about  the  z  axis  and  then  rotating  the  intermediate 
(x  ,  y\  zb  system  by  the  angle  e  about  the  x^axis. 

The  direction  cosines  between  the  spatial  and  cable  coordinate  systems  are 
obtained  as  the  elements  of  the  square  matrix  [a^  ,  given  by  the  product  of  the 
elementary  transformation  matrices  CO  and  £♦»] 
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These  direction  cosines  are  also  shown  in  Figure  1.  The  matrix  M  is  used 
to  express  the  components  of  a  vector  given  in  the  spatial  (x,  y,  z)  coordinate 
system  in  terms  of  the  cable  (X,  Y,  Z)  system,  while  the  matrix  0"  J  is  used 
to  express  the  components  of  a  vector  given  in  the  cable  (X,  Y,  Z)  systems  in 
terms  of  the  spatial  (x,  y,  z)  system.  Since  the  transformation  matrix  [A]  is 
orthogonal,  the  inverse  of  w  n  ,  is  simply  the  transpose  of  H  ■  [•'] 


[«"]  =  [«T] 


(4) 


CURRENT  PROFILE 

4 

As  in  Program  CAB3E  ,  the  current  profile  is  restricted  to  components 
in  the  horizontal  x  and  z  directions,  which  is  the  usual  case  for  ocean  current 
profiles.  However,  the  formulation  in  Program  CAB3E  is  generalized  to  allow 
the  current  magnitude  c  and  angle  with  the  x-axis,  oC  ,  to  have  steady-state  as 
well  as  time-varying  components,  as  follows 

~  (~o(y)  •+  ly) 

where  the  subscripts  o  and  1  denote  steady-state  and  time-varying  components, 

respectively,  f  is  the  frequency  of  the  current  variation 

c 

t  is  the  time 
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The  x  and  z  components  of  the  current,  c  and  c  ,  are  given  by 

x  " 

(7) 

(8) 

PRESCRIBED  SURFACE  MOTIONS 


Cv  -  C  OCX  o( 


Cj  C  c  5 in* 


As  in  the  case  of  the  Program  CABMOD,  the  program  allows  the  user  either 
to  prescribe  the  motion  at  the  surface,  as  in  a  towing  maneuver  case,  or  describe 
the  motion  by  means  of  differential  equations  for  a  surface  buoy. 

The  prescribed  motion  may  be  composed  of  a  series  of  sinusoidal  components 


in  (a)  the  x  and  y  directions,  or  (b)  the  x  and  z  directions. 

*  u*  (-irrht  +  <t>,  ) 


(a) 


i-l 

ti 


i-l  0 


(9a) 


(9b) 


(b) 


N  . 

x-  Z  t  +*i) 


L-l 

ti 


(10a) 


(10b) 


where  the  subscript  s  denotes  the  surface 

a  a  and  a  .  are  the  amplitudes  of  the  ith  sinusoidal  component 

xx  y  l  z  x 

In  the  x,  y,  and  z  directions,  respectively. 
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fj>  and  (j)^  are  the  frequency  and  phase  angle  of  the  ith  sinusoidal 
component,  respectively. 

Prescribed  motion  (a)  may  represent  a  surface  ship  in  waves  while  prescribed 
motion  (b)  may  represent  a  circular  towing  case. 

The  user  may  make  use  of  the  subroutine  PRESC  to  prescribe  more  general 
towing  maneuvers.  For  example,  a  build  up  in  towing  speed  V  from  0  to  If)  knots 
may  be  prescribed  according  to  the  formula 

V  -  10  X  1.6878  tanh  (t/200)  (11) 


SURFACE  BUOY  FORCES 

As  mentioned  In  the  Introduction,  surface  waves  and  rotational  modes  of  the 
buoy  are  considered  to  be  perturbations  to  the  overall  mean  motion  and  have  been 
neglected.  This  leads  to  a  number  of  considerable  simplifications.  First,  the 
exciting  forces  due  to  surface  waves  need  not  be  computed.  Secondly,  equations 
of  motion  need  not  be  written  for  the  rotational  modes,  thus  also  eliminating 
the  coupling  between  the  rotational  and  translational  modes.  Thirdly,  locations 
of  various  forces  need  not  be  computed  or  specified. 

For  the  sake  of  convenience,  it  is  further  assumed  that  the  buoy  is  axisym- 
metrlc  so  that  its  added  Inertia  and  drag  area  are  the  same  for  motion  in  any 
horizontal  direction.  This  Is  the  usual  case  for  the  surface  buoy  of  free- 
floating  cable  systems.  Under  the  above  assumptions,  the  resultant  forces  acting 
on  the  surface  buoy  are  as  follows: 
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w' ere  the  meanings  of  the  variables  are  as  follows:  I,  D,  G,  B,  WF,  T  * 
(inertia,  drag,  gravity,  buoyancy,  wind,  and  tension  forces,  respectively), 
~l,  ~]t  and~k  =  unit  vectors  in  the  x,  y,  and  z  directions,  respectively. 

>1  M  =  virtual  mass  in  the  horizontal  and  vertical  directions, 

SVH  ’  SVy 

respectively , 
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intermediate  hoik  forces 

The  forces  acting  on  the  bodies  placed  along  the  cable  are  somewhat  similar 
to  those  previously  derived  for  the  surface  buoy.  The  principal  differences 
occur  due  to  the  fact  that  the  bodies  have  a  constant  submerged  volume,  experience 
no  wind  force  and  have  cable  tension  forces  acting  above  and  below.  In  the  case 
of  the  lower  body,  a  special  provision  is  made  for  a  time-varying  thrust  which 
may  serve  to  alter  the  configuration  of  the  cable  system.  The  resultant  forces 
are  as  follows. 


T6 


X  *-  + 


f 'i ^ 
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. — ^ 

'  z^r^r^r  (^iAh +  (w«) 


<9,5  £  Co/i&^  £  2 
cp  As^ 


+  Men  -fK+ti  x  k  )  c>f*2L.  <0-5  <a  >2 

^  Cr»  Aftt* 


6-„  +  B„  =  W, 


6 1 


(19b) 


Lower  body  only: 

BT  -  6T 


caa  ( 2 iT-fr  t  t^r)?  +  *'•"  (2n^  ^i)  k  J 


where  the  subscript  B  denotes  the  body,  the  subscripts  a  and  b  refer  to  above 

and  below  the  body,  respectively,  W  is  the  weight  in  fluid  of  the  body,  BT, 

B 

f^,  and  <}>  are  the  magnitude,  frequency,  and  phase  of  the  time  varying  force 
acting  on  the  lower  body. 
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The  formulation  for  drag  force  given  in  Equation  (19a)  is  appropriate  for  a 
blunt  shape  such  as  a  cube  or  sphere.  The  formulation  given  in  Equation  (19b)  is 
more  appropriate  for  a  cylinder  or  disk  shape. 


CABLE  FORCES 

Following  Program  CABMOD,  the  continuous  cable  is  divided  into  a  number  of 
massless  straight  elastic  segments,  as  shown  in  Figure  2.  All  the  cable  forces 
are  equally  divided  between  the  two  nodes  at  the  ends  of  the  segment.  The  nodes 
may  represent  convenient  locations  to  represent  cable  shape  or  may  correspond  to 
locations  of  bodies  placed  along  the  cable,  including  the  surface  buoy  and  lower 
body. 

The  formulation  for  cable  forces  is  similar  to  the  formulation  in  Program 
CABMOD  with  three  principal  differences.  First,  a  somewhat  simpler  tension- 
strain  relation  is  used,  as  given  by  Equation  (2).  Secondly,  the  formulation 
is  extended  to  three  dimensions  using  the  coordinate  systems  shown  in  Figure  1. 
Thirdly,  added  inertia  terms  for  the  cable  are  neglected.  This  not  only  reduces 
the  number  of  terms  but  also  eliminates  coupling  of  the  equations  of  motion  for 
a  given  node. 

The  forces  are  given  for  a  cable  segment  of  lengths  .  Recall  that  these 
forces  must  be  equally  divided  between  the  nodes  at  the  ends  of  the  segment. 

The  two  simplest  forces  which  require  no  coordinate  transformation  are  the 

inertia  force  T  and  weight-in-fluid  force 

c  c 

Xt  =  Ml  yj  +  '}*■)  (23) 
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w c  =  wlH 


(24) 


where  the  subscript  c  denotes  the  cable 
JJi  Is  the  mass  per  unit  length 
AjJ  is  the  weight  per  unit  length  In  fluid 


The  tension  and  drag  forces  require  coordinate  transformations  and  will 
be  discussed  separately.  The  tension  force  T  is  given  in  Equation  (2)  as 


where  i  is  the  cable  length  at  T  *  T 

r  r 

t,H  is  the  change  in  length 

the  subscripts  %  and  u  refer  to  the  lower  and  upper  ends 
of  the  cable  segment,  respectively 


The  components  of  the  tension  in  the  spatial  x,  y,  and  z  directions  ate 
obtained  by  using  the  matrix  given  by  Equations  (3)  and  (4),  and  noting 

that  the  tension  force  for  the  case  where  the  cable  is  below  the  node,  ,  is 
directed  in  the  positive  Y  direction. 
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(26) 


which  leads  to 


T4 


-  oa  %  $in  ^  Tb 

On  ocn  Ty, 

5»”  %  ~t, 


(27a) 

(27b) 

(27c) 


cos  (y ,  Y)  -  ££  -  urx  9  <#*■  h,  < 

"  V 

cos  (z,  Y)  -  £1  -  5»r>&  ( 

H 

Equation  (?9a)  divided  by  Equation  (29b)  yields  the  following  equation  for  q>^ 

•  ’  L  %-•).)}  < 


(29b) 


(29c) 


while  Equation  (29c)  divided  by  Equation  (29b)  yields  the  equation  for  0 


The  normal  and  tangential  drag  forces  are  taken  to  be  respectively  pro¬ 
portional  to  the  squares  of  the  normal  and  tangential  components  of  the  fluid 
velocity  relative  to  the  cable  vr  given  by 

^  ^  — >  »  — ^  i  **♦  i  A 

Vr  ^  -Vc  +  c  -  0Cr  A.  t^rJ  4-  (? 


where 


Y  *  -^  + 

*  -~i 


c  and  c  are  respectively  the  x  and  z  components  of  current  given 

X  Z 

by  Equations  (5)  to  (8). 
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The  components  of  in  the  cable  X,  Y,  and  Z  directions  are  obtained  by  using  the 
matrix  [j^],  Equation  (3)  t  _ 


--  > 
Sr 


which  leads  to 


1TX  -  u*-4> +  si"^  > 


V s  ■  ** 


6  sin^#r  +  C»0«».^  A  +-  Sin© 


■\T  -  6in0 ir-  5/n0  un<j>vyr  +  cot-0  ^ 


(34a) 


(34b) 


(34c) 


The  drag  forces  D  ,  D  ,  and  0  in  the  X,  Y,  and  Z  directions  are  given  by 

LX  LY  L/C 


dcy  =  if  CT  H  ^  Ky  I 

D cl  z  ^  ^  ^  vr\  f  i rr\ 


(35a) 


(35b) 


(35c) 


where  d  is  the  cable  diameter 


C  is  the  normal  drag  coefficient 
D  6 

C  ,  is  the  tangential  drag  coefficient 
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In  order  :o  use  these  forces  in  the  equations  of  motion  it  is  necessary  to 
express  them  in  terms  of  the  spatial  directions  by  using  the  matrix 


--[A'']  « 


Da 


which  leads  to 


cast  4^  0 ^  -  6tr3.0sm4^  DcY45,ne  sin (37a) 


-  s)n  K  DCX  4coa'ec^<t  DCi^s^eux^\/D^ 


(37b) 


sin&  +  6w-  0  Pci  (37c) 


FINAL  EQUATIONS  OF  MOTION 

For  an  intermediate  node,  the  final  equations  of  motion  are  obtained  by 
using  the  body  forces  given  by  Equations  (18)  to  (21)  and  adding  to  them  the 
contributions  from  the  forces  of  the  cable  segments  above  and  below  the  node, 
Equations  (23)  to  (37),  resulting  in 

i6+i(ic4+icb)  =k J 

+  [Vi(Pc*+3a)]  +  [Zti]  ( 
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For  a  general  t  ^.e-dimensional  case,  the  above  vector  equation  leads  to  three 
uncoupled  scali -  second-order  differential  equations  for  the  x,  y,  and  z 
locations  of  the  node. 

For  the  lower  body,  the  quantities  Involving  the  subscript  b  are  equal  to 
zero  since  there  Is  no  cable  below  the  body.  In  addition,  there  is  a  thrust 
force  BT  acting  in  the  x  and  z  directions. 

For  the  surface  buoy,  the  vector  equation  of  motion  is  obtained  by  using  the 
buoy  forces  given  by  Equations  (12)  to  (17),  and  adding  to  them  the  contributions 
from  the  forces  of  the  cable  segment  below  the  buoy,  Equations  (23)  to  (37), 
resulting  in 

■ — »  (39) 

4  hi  Dcb)+  WF j  +  Tk 

DESCRIPTION  OF  COMPUTER  PROGRAM 

Trogram  CAB3DYN  consists  of  a  main  program  and  five  subroutines. 

MAIN  PROGRAM 

The  main  program  accepts  input  data  for  the  cable  system,  current  profile, 
initial  conditions,  and  data  for  either  a  surface  buoy  or  prescribed  motions 
according  to  Equations  (9)  or  (10).  Data  may  be  entered  In  either  English  or 
metric  units.  A  detailed  description  of  input  directions  is  given  in  the  next 
chapter. 

The  program  is  currently  written  to  accept  up  to  25  cable  segments  and  inter¬ 
mediate  bodies.  This  number  can  be  easily  increased  by  changing  a  few  DIMENSION 
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and  COMMON  statements,  but  it  should  be  noted  that  calculations  for  more  than  25 
nodes  are  likely  to  require  prohibitively  large  amounts  of  computer  time. 

The  main  program  prints  out  the  input  data,  and  then  calls  on  the  subroutine 
KUTMER  to  calculate  the  motions  of  the  system  at  prescribed  time  intervals.  The 
main  program  then  prints  out  the  calculated  displacements,  velocities,  accelera¬ 
tions,  and  tensions. 

SUBROUTINE  DYNA 

This  subroutine  defines  the  vector  differential  equation  of  motion  (38)  for 
each  of  the  N  nodes.  The  subroutine  also  defines  the  vector  differential  equation 
of  motion  (39)  for  the  case  of  a  surface  buoy,  or  the  prescribed  motion  according 
to  Equations  (9)  or  (10).  In  the  case  of  a  surface  buoy,  this  subroutine  calls 
the  subroutine  XLINEAR  to  furnish  the  values  of  buoy  volume  and  drag  area  as  a 
function  of  buoy  submergence.  In  the  case  of  a  prescribed  motion  not  covered  by 
Equations  (9)  or  (10),  the  subroutine  calls  the  subroutine  PRESC  for  the  motions. 

SUBROUTINE  CUR 

This  subroutine  furnishes  the  x  and  z  components  of  the  time-varying  current 
profile.  For  a  given  value  of  the  vertical  distance  y,  the  subroutine  linearly 
interpolate  between  the  quantities  c  ,  ot  ,  c.  ,  and  ex,  ,  defined  in  Equations 
(5)  and  (6),  which  are  input  as  a  function  of  y.  The  subroutine  then  calculates 
the  x  and  z  components  of  current  according  to  Equations  (7)  and  (8).  For  cases 
where  the  given  value  of  y  is  greater  (less)  than  the  largest  (smallest)  value  of 
y  which  is  read  in,  the  subroutine  takes  the  values  of  ce  ,0/^,  c(  ,  and  ,  to 
be  the  values  at  the  largest  (smallest)  value  of  y  which  is  read  in. 
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SUBROUTINE  KUTMER 


Thi3  subroutine  uses  the  Kutta-Merson  method  to  numerically  Integrate  the 
differential  equations  of  motion  defined  in  subroutine  DYNA.  The  subroutine 
automatically  reduces  the  integration  step  size  uni  11  i'v  specified  error 
criteron  of  0.0001  is  satisfied. 

SUBROUTINE  XLINEAR 

This  subroutine  furnishes  the  volume  V  and  drag  area  C^A^  as  a  function 
of  buoy  submergence  H.  The  subroutine  linearly  interpolates  between  the  input 
values  of  volume  and  drag  area  which  are  read  in  as  a  function  of  submergence. 

For  cases  where  the  given  submergence  is  greater  (less)  than  the  largest 
(smallest)  value  of  submergence  whLoh  is  read  in,  the  subroutine  takes  the  volume 
and  drag  area  to  be  the  value  at  the  largest  (smallest)  value  of  submergence 
which  is  read  in. 

SUBROUTINE  PRESC 

The  user  may  use  this  subroutine  to  prescribe  an  arbitrary  motion  which  is 
outside  the  form  covered  by  Equations  (9)  and  (10). 

PROGRAM  STORAGE  AND  TIME  REQUIREMENTS 

On  the  CDC  6700  computer  currently  in  use  at  the  Center,  Program  CAB3DYN 
requires  a  memory  of  approximately  41,000  octal  words  to  load,  30,400  octal  words 
to  execute,  and  approximately  18  seconds  to  compile.  By  comparison,  Program 
CABM0D  requires  a  memory  of  approximately  52,500  octal  words  to  load,  41,000 
octal  words  to  execute,  and  approximately  30  seconds  to  compile.  Thus,  Program 
CAB3DYN  is  a  simple  program,  despite  its  three-dimensional  capability.  This 
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is  largely  due  to  the  previously  mentioned  neglect  of  surface  waves,  cable 
added  inertia,  and  buoy  rotational  modes. 


As  in  the  case  of  Program  CABMOD,  execution  time  depends  largely  on  the 

value  of  the  elastic  constant  C  ,  the  number  of  cable  nodes,  and  the  amount 

1 

amount  of  time  over  which  dynamic  motions  are  desired.  For  the  sample  problems 
given  in  this  report,  execution  time  for  Program  CAB3DYN  ranges  from  1.2  to  2.5 
times  the  corresponding  time  for  Program  CABMOD.  In  these  problems,  the  ratio  of 
execution  time  to  real  time,  ET/RT,  ranges  from  0.6  to  3.3. 


INPUT  INSTRUCTIONS 


READ  STATEMENTS 


Input  data  are  entered  into  the  program  by  means  of  the  following  READ 
statements  in  the  MAIN  program.  Wherever  possible,  the  names  and  definitions 


of  Input  variables  are  matched  to  those  contained  in  Program  CABMOD. 

READ(3, 1 )  HCASES 
BO  1000  HOI'HCASES 
READ ( 5 , 301 )  TITLE 
301  F0RHAT  <  20A4 ) 

READ<5,1>  NSMCAB,NCUMTRC,ISK 
RE AD < 5,2)  (FSH(K) ,K=1 ,NSM) 

2  FORMAT ( 8F 10.4) 

3  F0RHAT<8F10.2) 

4  FORMAT  <  8F 1 0 . 6 ) 

5  FORMAT  <  8F 1 0 . 0  > 

READ<5,2>  (AXSN(K) , K = 1 ,RSH) 

READ < 5 1 2 )  ( AYSM < K ) , K  =  1 ,NSM) 

READ < 5,2)  (AZSH(K) ,Ks1 ,NSH) 

READ (5,2)  <FI0SH(K) ,K=1 ,NSN) 

READ (5,2)  RH0 fTMIMyCCFR 

READ (5,2)  BTM , BTF , BTPD 

READ ( 5 f 2 )  TINV1 ,DT1 ,T0TT,DT2,TBYHX 

RE AD < 5 , 3 )  (FLC(K) ,K*1 , NC AB ) 

READ ( 3,2)  (DCI (K) ,K=1 , NC AB ) 

READ( 5, 2)  (CDN(K) ,K=1 , NC AB ) 

READ (5,2)  (CDT(K),K=1 , NCAB ) 

RE AD < 5,2)  < UC < K ) ,K=1 .NCAB) 
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READ < 5 p 4 )  (CH(K)  ,K=1 , NCAB) 

RE  AD ( 5 , 3 )  (TREF(K) ,K=1 , NC AB ) 

READ < 5 , 5 )  (Cl (K) ,K=1 r NCAB ) 

READ < 5 r 2 )  <UBD(K) , K-1 f NCAB) 

READ (5,2)  (CDABX(K),K=1 , NCAB ) 

READ <  5 f 2  >  < CDABY < K ) , K= 1 f NCAB ) 

READ < 5 f 2 >  (XHBV(K) ,K=1 ,NCA8) 

READ (5,2)  ( YMBV < K ) » K  =  1 f NCAB ) 

READ(5f 3)  (YY ( I) , 1=1 ,NCUR) 

RE AD < 5 r 3 )  (CCZK(K) ,K*1 , NCUR ) 

READ ( S , 3 )  (CC1K(K),K=1,NCUR) 

READ (5,3)  (CCZAD(K) ,K*1 f NCUR ) 

READ (5,3)  (CC1AD(K),K*1 ,NCUR) 

READ < 5 1 2 )  (PHID(I), 1  =  1, NCAB) 

READ (5,2)  (THETAD(I), 1=1, NCAB) 

READ ( S , 3 )  (TENI(I), 1  =  1, NCAB) 

READ (5,2)  (XPI(l), 1*1, NCAB) 

READ (5,2)  (YPI(I), 1=1, NCAB) 

READ ( S , 2 )  (ZPI(I),  1*1, NCAB) 

IF  (FSH( 1 ) .LE.2000. )  GO  TO  220 
READ (5,2)  XSI ,YSI ,ZSI ,XPSI, YPSI ,ZPSI ,UFN,UFAC 
READ ( S y 2 )  CDASBY,UASB,XHSBV, YNSBV,RTX,RTY,RTZ,DFTLIH 
220  CONTINUE 


The  corresponding  FORMAT  statements  are: 

1  FORMAT  (2413) 

2  FORMAT  (8F10.4) 

3  FORMAT  (8F10.2) 

4  FORMAT  (8F10.6) 

5  FORMAT  (8F10.0) 

301  FORMAT  (20A4) 


DEFINITION  OF  INPUT  VARIABLES 


NCASES 

Number 

of 

cases,  NCASES  >  1 

TITLE 

Title 

NSM1  * 

Number 

of 

surface  motion  components 

NCAB 

Number 

of 

cable  segments,  2  <  NCAB 

1_<  NSM 
<  25 


< 
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*  Superscripts  refer  to  explanatory  notes  which  begin  on  page  27 
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NCUR 

MTRC 

ISM 


FSM(K)1 

AXSM(K) 

AYSM(K) 1 . 
AZSM(K) 

FIDSM(K)1 

RHO 

TMIN 

CCFR 

BTM,  BTF, 
BTPD 

TINV1 

DTI 

TOTT 

DT2 


Number  of  current  profile  points,  2<_  NCUR  <_  10 
MTRC  0  if  input  data  are  enetered  in  English  units; 

MTRC  1  if  input  data  are  entered  in  metric  units** 

ISm£0  if  surface  motions  are  prescribed  in  x  and  y  directions, 
Equation  (Q);  ISM  ^ 1  if  surface  motions  are  prescribed  in  x  and 
z  directions,  Equation  (10) 


Fluid  <&ns*ity  in  slugs /feet^or  kilograms  ^ meters^ 


7SM 


Minimum  tension  which  can  be  supported  by  cable  in  pounds 
or  Newtons;  TMXN  *  0  for  a  flexible  cable 
Frequency  of  current  variation  in  cycles  per  second 
Bottom  thrust  BT  ■  BTM# 


•  *  ^ 

+■  *»  n  (lift  &7F  *t  -t  )  k  j 


Initial  time  interval  in  seconds  for  motion  calculations 
Time  step  in  seconds  for  which  printout  is  desired 
0  <_  t  _<  TINV1 

Total  time  in  seconds  for  which  motion  calculations  are  desired 
Time  step  in  seconds  for  which  printout  is  desired  for 
TINV1  <  t  <  TOTT 


**  All  data  must  be  input  in  consistent  units. 
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TBYMX 

FLC(K) 

DCI(K) 

CDN(K) 

CDT(K) 

WC(K) 

CM(K) 

TREF(K) 

C1(K) 

WBD(K) 

CDABX(K) , 
CDABY(K) 

XMBV(K) 

YMBV(K) 

YY(I) 

CCZK(K) 

CC1K(K) 

CCZAO(K) 

CCIAD(K) 

PHID(I) 

THETAD(I) 

TENI(I) 

XPI(I),  YPI(I) 
ZPI(I) 


Maximum  absolute  value  of  tension  in  cable  just  below  buoy 
in  pounds  or  Newtons 

Length  of  Kth  cable  segment  in  feet  or  meters 
Diameter  of  Kth  cable  segment  in  inches  or  centimeters 
Normal  drag  coefficient  of  Kth  cable  segment 
Tangential  drag  coefficient  of  Kth  cable  segment 
Weight  in  fluid  of  Kth  cable  segment  in  pounds/foot  or 
Newtons/meter 

Mass  of  Kth  cable  segment  in  slugs/foot  or  kilograras/meter 

Reference  tension  of  Kth  cable  segment  in  pounds  or  Newtons 

Elastic  constant  of  Kth  cable  segment  in  pounds  or  Newtons 

Weight  in  fluid  of  Kth  body  in  pounds  or  Newtons 

Drag  area  of  Kth  body  in  feet?  or  meters  ^  for 
flow  in  (horizontal, vertical)  directions 

Virtual  mass  (mass  t  added  mass)  of  Kth  body 

in  slugs  or  kilograms  for  motion  in  (horizontal, vertical) 

directions 

Value  of  y  in  feet  or  meters 

Magnitudes  ( )  of  (steady,  time-varying)  component  of 
current  in  knots  or  meters/second  at  y  *  YY(K),  see  Equations 
(5)  to  (8) 

\ngles  (ofC/ )  measured  from  x-axis  of  (steady,  time-varying) 
component  of  current  in  degrees  at  y  *  YY(R),  see  Equations  (5) 
to  (8) 

Initial  value  of  $  of  Ith  cable  segment  in  degrees 

Initial  value  of  8  of  Ith  cable  segment  in  degrees 

Initial  value  of  tension  of  Ith  cable  segment  in  pounds  or  Newtons 

Initial  values  of  (x,  y,  z)  of  Ith  node  in  feet/second  or 
meters/second 
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DEFINITION  OF  INPUT  VARIABLES  FOR  SURFACE  BUOY,  (FSM(l)  >  2000) 


XSI(I),  YSI(I),  Initial  values  of  (x,  y,  z)  of  buoy  center  in  feet  or  meters 
ZSI(I) 


XPSI(I),  YPSI(I),  Initial  values  of  (x,  y,  z)  of  buoy  center  in  feet/secrml 


or 


ZPSI(I) 

WFM 

WFAD 

CDASBY 

WASB 

XMSBV 

YMSBV 

RTX.RTY, 

RTZ 

DFTLIM 


meters/secnn 1 

Magnitude  of  wind  force  in  pounds 

Angle  with  x-axis  of  wind  force  in  degrees 

2  2 

Drag  area  for  y-  direction  in  feet  or  meters 
Weight  in  air  in  pounds  or  Newtons 

Virtual  mass  (mass  +  added  mass)  for  motion  in  (horizontal, 
vertical)  directions  in  slugs  or  kilograms 

Distance  of  cable  attachment  point  from  buoy  center  in  (x,  y,  z) 
directions  in  feet  or  meters 

Limiting  value  of  draft  in  feet  or  meters  for  which  the  pro¬ 
gram  calculates  the  percent  of  time  that  the  draft  exceeds  this 
value . 


EXPLANATORY  NOTES 

1.  For  FSM(l)  >  2000.,  the  program  accepts  input  data  for  a  surface  buoy. 

3  3 

In  this  case,  AXSM(K)  is  the  submerged  volume  of  the  buoy  in  feet  or  meters 
at  vertical  distance  AYSM(K)  in  feet  or  meters  from  the  buoy  center,  defined  as 
the  origin  of  the  local  buoy  coordinate  system.  Input  data  is  entered  starting 
with  the  bottom  of  the  buoy  and  progressively  moving  up  so  that  AYSM(l)  is  the 
largest  algebraic  value  and  AYSM(NSM)  is  the  smallest  algebraic  value.  Also, 

AYSM(l)  -  AVS'i('JS'i)  =*  total  buoy  length  1^  .  Also,  AZSM(K)  is  AgH,  the  submerged 

2  2 

drag  area  in  the  horizontal  direction  in  feet  or  meters  at  vertical  distance 
AYSM(K)  from  the  buoy  center.  The  FiDSM(K)  may  take  on  any  values,  such  as, 
say,  0. 
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For  FSM(l)  <0.,  the  program  calls  the  subroutine  PRESC  for  the  pre¬ 
scribed  surface  motions.  This  subroutine  is  used  by  the  user  to  program  pre¬ 
scribed  motions  which  are  not  covered  by  Equations  (9)  and  (10). 

2.  For  FSM(l)  >  2000.,  the  program  calculates  the  draft  11  ,  required  to 
support  the  weight  in  air  of  the  buoy,  and  the  weight  in  fluid  of  the  lower 
cables  and  bodies.  The  displacement  YYSI  of  the  buoy  center  from  the  ocean 
surface  for  this  draft  is  given  by 

YYSI  =  Hs-AY5M  (i)  (40) 

O 

For  an  input  value  of  YSI  <  1.0  X  1(T  ,  the  program  sets  the  initial  value 

of  YSI  “  YYSI.  This  feature  is  specially  useful  for  those  cases  where  the  cable 
starts  from  an  initially  vertical  position. 

SAMPLE  PROBLEMS 

Two  sample  problems  are  given  to  illustrate  usage  and  results  calculated 
by  the  program. 

FREE-FLOATING  BUOY-CABLE  PROBLEM 
Steady  Current 

The  first  problem,  shown  in  Figure  3,  involves  the  deployment  behavior  of  an 
initially  vertical  100-foot  (30.5  m)  free-floating  buoy-cable  system,  shown  in 
Figure  3,  in  the  presence  of  a  current  profile  whose  magnitude  varies  linearly 
from  2  knots  at  the  surface  to  zero  at  y  =  100  feet  (30.5  m) ,  as  shown  in  Figure  3. 
Calculations  were  made  for  a  series  of  variations  of  current  direction  o(  .  Table  1 
shows  for  each  direction  variation,  the  magnitude  V^and  angle  o^of  the  final 


steady-state  drift  velocity,  the  depth  y  and  horizontal  excursions  x  and  z 

L  L  L 

of  the  lower  body  measured  from  the  surface  buoy,  the  real  time  t  required 
to  reach  the  final  steady-state  configuration,  and  the  ratio  of  execution 
time  to  real  time,  ET/RT.  The  table  also  shows,  for  comparison  purposes,  two 
sets  of  steady-state  results  calculated  by  the  two-dimensional  Program  CABMOD. 
The  first  set  is  obtained  by  directly  integrating  the  differential  equations 
of  equilibrium  and  iterating  the  drift  velocity  and  buoy  draft  until  f he  cable 
system  is  in  equilibrium.  This  approach  will  yield  more  accurate  steady-state 
results  than  the  approach  used  in  Program  CAB3DYN  since  the  cable  is  taken  to 
be  •'ontinuous  rather  than  segmented.  The  second  set  of  results  calculated  by 
Program  CABMOD  is  obtained  in  a  manner  similar  to  that  of  Program  CAB3DYN, 
where  the  cable  is  segmented  and  the  deployment  motion  is  calculated  from  an 
initially  vertical  position  to  the  final  steady-state  shape. 

The  results  shown  In  the  table  may  be  used  to  check  the  accuracy  of  the 
formulation  and  program  logic  in  Program  CAB3DYN  in  three  ways.  First,  Cases 
1  and  2  represent  :he  same  unidirectional  current  profile  with  reference 
direction  shifted  by  45  degrees.  The  table  shows  that  the  only  difference  in 
the  calculated  results  for  the  two  cases  is  a  shift  of  45  degrees  in  the  drift 
velocity  and  horizontal  displacements  x  and  z  .  This  shows  that  the  cal- 

Li  L 

culation  procedure  for  the  horizontal  x  and  z  displacements  are  equivalent. 
Secondly,  a  comparison  of  Cases  1  and  6  shows  that  Program  CAB3DYN  predicts 
values  of  V  and  x  which  respectively  differ  by  2.7  and  0.4  percent  from  the 

L)  Li 

exact  steady-state  formulation  represented  by  Case  6.  This  comparison  indicates 
the  type  of  error  Incurred  when  the  continuous  cable  is  represented  by  four 
segments.  Perhaps  the  most  direct  check  is  to  compare  Cases  1  and  7,  since 
they  represent  the  final  steady-state  results  calculated  by  Programs  CAB3DYN  an< 
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CABMOD,  respectively,  starting  from  identical  initial  states  and  identical 
segmented  representations  of  the  cable.  Table  1  shows  that  the  agreement  is 
only  to  two  significant  places.  Table  1  also  shows  that  the  CABMOD  results 
for  Vp  and  x^  respectively  show  oscillations  of  +  0.01  ft/sec  (+_  0.003  m/sec) 
and  +_  0.1  ft  (+  0.03  m)  about  the  average  values.  The  reason  for  both  of  the 
above  phenomena  is  that  Program  CABMOD  models  buoy  pitch  which  still  has  not 
damped  out  at  t  *  40  sec,  while  Program  CAB3DYN  neglects  all  rotational  modes 
of  the  buoy.  The  next  problem  gives  a  more  valid  comparison  since  the  surface 
motion  is  identically  modeled  by  both  programs.  The  results  do  show  that  the 
presence  of  buoy  pitch  has  relatively  little  effect  on  the  average  results. 

Cases  3,  4,  and  5  show  results  for  three-dimensional  configurations  where  the 
current  changes  in  direction  with  depth.  It  is  of  interest  to  note  that  the  effect 
of  change  in  direction  Is  to  slow  drift  velocity  by  as  much  as  15%  from  the  two- 
dimensional  value  given  by  Cases  1  and  2.  Furthermore,  the  resultant  horizontal 
displacement  K  +  >L  for  the  three-dimensional  Case  5  is  20  percent  larger 
than  the  two-dimensional  value.  Thus,  the  use  of  a  unidirectional  current  need 
not  predict  the  largest  horizontal  displacement. 

The  table  shows  that  the  time  required  to  deploy  from  an  initially  vertical 
position  to  the  final  steady-state  configuration  is  approximately  40  seconds  for 
all  cases.  The  ratio  of  execution  time  to  real  time  is  approximately  3.2  for 
Program  CAB3DYN  and  1.2  for  Program  CABMOD. 

It  is  of  interest  to  know  the  effect  of  cable  added  inertia,  which  is  neg¬ 
lected  in  Program  CAB3DYN,  on  the  cable  motions  during  deployment.  The  motions 
predicted  by  Programs  CABMOD  and  CAB3DYN  are  not  directly  comparable  due  to 
to  their  previously  noted  difference  in  the  modeling  of  the  surface  buoy. 
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However,  it  should  be  noted  that  Program  CABMOD  can  be  run  with  and  without 
cable  added  inertia  by  setting  the  input  variable  AMC  equal  to  1.  and  0., 
respectively.  Table  2  shows  values  of  the  velocities  x  at  the  cable  nodes 
computed  by  Program  CABMOD  for  AMC  *  1  and  AMC  =*  0  for  different  times  t  during 
deployment.  In  the  table,  Node  0  represents  the  point  of  attachment  of  the 
cable  to  the  buoy  while  Node  4  represents  the  lower  weight.  The  table  also 
shows  the  ratio  ET/RT  for  AMC  *  1  and  AMC  *  0.  The  results  show  that  cable 
added  inertia  has  little  effect  on  the  motions  since  the  maximum  difference 
between  the  results  for  AMC  =  1  and  AMC  =»  0  is  only  0.004  ft/sec  (0.0012  m/sec). 

As  expected,  the  largest  differences  occur  at  the  beginning  of  the  deployment 
when  inertia  effects  are  the  largest.  The  table  also  shows  that  ET/RT  is  1.2 
for  the  AMC  *  1  case  compared  to  1.1  for  the  AMC  =  0  case. 

Comments  on  Efficient  Use  of  Program  for  Slowly  Varying  Current 

The  results  contained  in  Table  1  offer  suggestions  on  efficient  ways  of  using 
Program  CAB3DYN  for  the  following  problem  which  otherwise  would  require  prohibitive 
amounts  of  computer  time.  Suppose  it  is  of  interest  to  solve  the  above  problem 
with  the  following  two  differences.  First,  the  current  is  no  longer  steady  but 
exhibits  a  sinusoidal  variation  with  time  with  a  period  of,  say,  three  hours. 
Secondly,  the  shape  of  the  cable  is  desired  at  20  nodes  instead  of  only  4  nodes. 

A  straightforward  use  of  the  program  would  model  the  cable  with  20  nodes  and 
compute  the  motion  of  the  cable  over  three  hours.  Table  1  indicates  that 
execution  time  is  approximately  3.2  times  real  time  for  4  nodes.  Reference  1 
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Indicates  that  execution  time  is  proportional  to  the  square  cf  the  number  of 
nodes  NCAB 


£T  J  lH CAB?  (4l) 

Thus,  the  execution  time  ET  required  to  calculate  motions  for  three  hours 
for  20  nodes  would  be  approximately 

ET  ~  3. Z  X  (\)  X  3  ^  2  40  hours 

J  (42) 

Assuming  that  such  an  amount  of  computer  time  were  available,  the  cost  using 
the  lowest  priority  available  at  the  Center,  block  time  priority  PO,  would 
be  approximately 


Cost  =  240  X  $125  -  $  30,000 


(43) 


However,  the  results  shown  in  Table  1  suggest  alternate  ways  of  obtaining 
the  desired  results  at  much  lower  computer  cost.  Table  1  shows  that  the  cable 
requires  only  40  seconds  to  deploy  from  an  initial  vertical  state  to  the  steady- 
state  configuration  for  a  given  current  profile.  Since  this  is  a  much  smaller 
time  scale  than  the  period  of  the  current  variation,  the  cable  system  at  a  given 
time  is  essentially  in  equilibrium  with  the  instantaneous  current  profile. 

The  table  also  shows  that  the  steady-state  drift  velocity  obtained  by  using  only 
four  nodes  is  only  a  few  percent  different  from  the  exact  value  obtained  oy  taking 
the  cable  to  be  continuous.  These  facts  suggest  the  following  alternate  approach. 


Use  Program  CAB3DYN  with  relatively  few  cable  segments  to  calculate  the  steady- 
state  drift  velocity  at,  say,  20  equally  spaced  values  of  time  t  covering 
the  current  profile  cycle.  If  the  cable  is  started  from  a  vertical  position 
each  time,  deployment  time  would  be  approximately  40  seconds.  On  the  other  hand, 
if  the  more  efficient  method  of  using  the  steady-state  values  from  the  preceding 
time  as  the  initial  conditions  are  adopted,  deployment  time  to  reach  the  new 
steady-state  configuration  will  be  substantially  reduced  to  perhaps  a  few  seconds. 
The  20  calculated  drift  velocity  vectors  $  may  then  be  used  to  generate  20 
resultant  current  profiles  c^  given  by 

cr^J  -  (44) 

These  current  profiles  along  with  the  calculated  values  of  buoy  draft  may  then  be 
input  into  the  thr«e-dimensional  steady-state  Program  CAB3E^  ,  which  calculates 
the  configuration  at  any  number  of  positions  along  the  cable  at  the  cost  of  only 
a  few  seconds  per  case.  Thus,  even  if  one  assumes  that  deployment  time  is  40 
seconds  for  each  current  profile,  total  execution  time  ET  for  this  alternate 
approach  is  given  by 


ET  z  20  x  40  -  800  sec  <45) 

The  cost  using  overnight  priority  P2  would  be  approximately 

Cost  =■  800  X  $0.06  =  $48  (46) 

Thus,  the  computer  time  and  cost  are  substantially  lower  than  the  straight¬ 
forward  approach  while  still  yielding  results  of  comparable  accuracy. 
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TOWING  CABLE  PROBLEM 


The  second  problem,  shown  in  Figure  4,  involves  the  deployment  behavior 
of  an  initially  vertical  500-  ft  (152.4  m)  cable  towed  at  various  speeds. 
Calculations  were  made  for  the  cable  with  and  without  the  neutrally  buoyant 
array  shown  in  Figure  4.  The  cable  is  modeled  by  5  equal  segments. 

In  addition  to  Programs  CAB3DYN  and  CABMOD,  Program  CABTRAN  was  also 
used  to  calculate  the  deployment  motion.  In  this  program,  the  lower  array 
is  essentially  lumped  with  the  last  cable  segment.  In  order  to  make  the 
results  calculated  by  Programs  CAB3DYN  and  CABMOD  equivalent  to  those  calculated 
by  Program  CABTRAN,  the  array  drag  was  accounted  for  by  increasing  the  diameter 
and  drag  coefficient  C^  of  the  last  segment. 

For  cases  of  the  cable  towing  the  neutrally  buoyant  array,  Programs  CABMOD 
and  CAB3DYN  predict  converged  steady-state  results  for  only  the  1-knot  case. 

Table  3  shows  four  setB  of  calculated  result  s  for  th.-  <  ijle  of  each  cable  seg¬ 
ment  at  various  times  t  during  the  deployment  for  this  case.  The  table  shows 
results  calculated  by  Program  CABTRAN,  Program  CAB3DYN,  and  Program  CABMOD  with 
and  without  added  cable  inertia  effects.  At  t  *  1,500  sec,  when  the  deploynent 
is  complete,  the  table  also  shows  the  steady-state  cable  angle  calculated  by 
'’rogram  CABMOD  taking  the  cable  to  be  continuous.  Finally,  the  table  also  shows 
the  ratio  of  execution  time  to  real  time  for  each  of  the  programs. 

The  table  shows  that  the  cable  angles  calculated  by  Program  CABMOD,  with 
md  without  added  cable  inertia,  and  by  Program  CAB3DYN  all  agree  to  within  0.03 
degrees  during  the  entire  deployment  motion.  As  in  the  towing  cable  problem,  the 
results  show  that  cable  added  inertia  effects,  which  are  neglected  in  Program 
CAB3DYN,  have  little  effect  on  the  motion.  The  results  also  show  that  the  formu- 
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lation  and  logic  In  Programs  CABMOD  and  CAB3DYN  for  two-dimensional  cases  are 
essentially  equivalent.  The  results  for  Problem  1  have  already  shown  that 
the  formulations  in  Program  CAB3DYN  for  the  x  and  z  directions  are  equivalent. 

The  table  shows  that  the  cable  angles  predicted  by  Program  CABTRAN  for  the 
upper  segments  agree  with  the  CABMOD  and  CAB3DYN  results  to  within  two  degrees 
during  the  entire  deployment.  However,  there  are  large  differences  for  the  two 
lower  segments.  In  particular,  there  are  differences  of  ft  and  17  degrees  at 
t  ■  1,500  seconds.  The  table  shows  that  the  CABMOD  and  CAB3DYN  angles  at  t  • 

1,500  seconds  agree  to  within  three  degrees  with  the  exact  steady-statu  ivsults, 
where  the  cable  is  taken  to  be  continuous.  The  principal  error  in  the  configur¬ 
ation  calculated  by  Program  CABTRAN  is  that  there  is  too  large  a  jump  in  the 
angle  of  the  two  lower  segments  which  is  opposite  to  the  trend  along  the  rest  of 
the  cable. 

In  view  of  this  incorrect  steady-state  behavior,  Program  CABTRAN  was  exer¬ 
cised  for  a  series  of  towing  speeds  and  number  of  cable  segments  NCAB.  Table  4 
shows  the  steady-state  cable  angles,  obtained  at  t  *  1,500  seconds  for  towing 
speeds  =  1,  3,  5  and  10  knots  and  values  of  NCAB  ■*  5,  10  and  20  segments. 

The  table  also  lists  the  ratio  ET/RT  for  each  value  of  NCAB. 

The  table  shows  that  for  a  given  speed,  the  type  of  error  does  not  change  with 

number  of  segments.  For  the  case  of  1  knot,  the  previously  observed  large  Jump  in 

angle  between  the  two  lower  segments,  which  is  opposite  to  the  trend  along  the  rest 

of  the  cable,  occurs  again  for  NCAB  •  10  and  20.  For  V  >  3  knots,  the  type  of 

S  — 

error  Is  somewhat  different.  For  all  three  values  of  NCAB,  the  cable  angle  oscil¬ 
lates  along  the  cable,  with  the  largest  oscillation  occurring  in  the  lower  seg¬ 
ments.  The  correct  steady-state  configuration  would  exhibit  a  monotonic  variation 
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of  angle  along  the  cable. 

For  V  >  3  knots,  Programs  CAB3DYN  and  CABMOD  are  not  able  to  predict 

O 

converged  steady-state  solutions.  The  programs  predict  large  oscillations  in  the 
inclination  of  the  lowest  segment,  even  at  t  «  2,000  seconds.  The  problem  appears 
to  be  the  shallow  angle  of  the  lowest  segment,  which  contains  the  drag  of  the 
neutrally  buoyant  array.  Computer  runs  for  cases  of  no  array  and  a  heavy  array 
do  not  exhibit  the  above  instability  behavior. 

Due  to  the  difficulties  encountered  by  all  three  programs  for  the  neutrally 
buoyant  array,  runs  were  made  using  Program  CABMOD,  without  added  cable  inertia, 
and  Program  CABTRAN  for  the  case  of  no  array  for  several  towing  speeds  . 

Tables  5  and  6  show  the  cable  angles  calculated  by  both  programs  at  various  times 

during  deployment  for  V  =  1  and  5  knots,  respectively.  Program  CABMOD  was  used 

b 

instead  of  Program  CAB3DYN  due  to  its  somewhat  lower  execution  time,  see  Table 
3.  As  shown  in  Table  3,  the  cable  angles  calculated  by  Program  CABMOD  agree  to 
within  0.02  degree  with  corresponding  values  calculated  by  Program  CAB3DYN. 

Table  5  shows  that  for  »  1  knot,  the  calculated  cable  angles  for  all  the 
segments  agree  to  within  0.2  deg  during  the  entire  deployment.  Table  6  shows  that 
for  Vg  =*  5  knots,  the  agreement  is  somewhat  worse,  with  differences  which  are 
typically  one  to  two  degrees.  Table  5  also  shows  that  both  programs  predict  some 
overshooting  past  the  steady-state  angles.  For  both  towing  speeds,  the  calculated 
steady-state  angles,  which  correspond  to  a  straight  line,  agree  to  within  0.01  deg. 


SUMMARY 

On  the  CDC  6700  computer,  Program  CAB3DYN  requires  approximately  41,000 
octal  words  to  load,  30,400  octal  words  to  execute,  and  18  seconds  to  compile. 
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For  two  sample  problems  Involving  the  deployment  of  a  free-floating  buoy-cable 
and  a  towing  cable,  execution  time  ranges  from  1.2  to  2.5  times  the  corresponding 
time  required  by  the  two-dimensional  Program  CABMOD.  The  ratio  of  execution  time 
to  real  time  ranges  from  0.6  to  3.3. 

For  both  problems,  the  results  show  that  cable  added  inertia,  which  is  neg¬ 
lected  in  Program  CAB3DYN,  has  little  effect  on  the  motions  during  deployment. 

The  accuracy  of  the  formulation  in  Program  CAB3DYN  is  verified  by  noting  the 
agreement  with  Program  CABMOD  for  a  towing  cable  case,  and  the  equivalence  of  the 
results  calculated  by  Program  CAB3DYN  for  two  unidirectional  currents. 

The  results  calculated  by  Program  CAB3DYN  for  the  free-floating  buoy-cable 
problem  agree  only  in  an  average  sense  with  those  calculated  by  Program  CABMOD. 
This  is  due  to  the  presence  of  buoy  pitch  terms  contained  in  Program  CABMOD, 
which  are  neglected  in  Program  CAB3DYN.  It  is  shown  that  the  steady-state 
results  calculated  by  modeling  the  deployment  motion  with  four  cable  segments 
differs  by  only  a  few  percent  from  exact  results  obtained  by  considering  the 
cable  to  be  continuous.  An  efficient  method  is  indicated  for  using  Program 
CAB3DYN,  in  conjunction  with  the  three-dimensional  steady-state  Program  CAB3E , 
to  calculate  the  motions  and  detailed  shape  of  a  free-floating  cable  system  in 
the  presence  of  a  slowly  varying  current.  The  method  essentially  consists  of 
using  Program  CAB3DYN  with  relatively  few  segments  to  calculate  the  steady- 
state  draft  of  the  surface  buoy  and  current  profile  relative  to  the  cable  system 
at  a  series  of  given  times  during  the  current  cycle.  The  draft  and  relative 
current  profile  are  then  read  Into  Program  CAB3E  to  calculate  the  cable  config¬ 
uration  at  any  number  of  positions  along  the  cable. 
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For  cases  where  the  cable  is  towing  a  neutrally  buoyant  array,  Programs 
CABMOD  and  CAB3DYN  predict  unstable  motions  of  the  lower  segment  at  higher  towing 
speeds.  This  instability  is  apparently  due  to  the  shallow  angle  of  the  lower 
segment  since  it  disappears  for  cases  of  no  array  and  a  heavy  array.  Program 
CABTRAN,  which  neglects  all  Inertia  effects,  was  also  used  to  calculate  results 
for  the  towing  cable  problem.  This  program  predicts  incorrect  steady-state 
configurations  for  the  case  of  the  neutrally  buoyant  lower  array.  For  cases  of 
cable  alone,  the  cable  angles  predicted  by  Programs  CABMOD  and  CAB3DYN  usually 
agree  to  within  one  or  two  degrees  with  those  predicted  by  Program  CABTRAN 
during  the  entire  deployment  cycle. 
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Ocean  Surface 


cos  (X,x)  =  cos  <£y 
cos  (Y,x)  =  -cos  0  singly 
cos  (Z,x)  =  sin  9  sin  <£y 


cos  (X,  y)  -  sin  <£y 
cos  (Y ,y)  =  cos  8  cos  <£y 
cos  (Z,y)  =  -sin  0  cos^y 


cos  (X,z)  =  0 
cos  (Y,z)  =  sin  8 
cos  (Z,z)  =  cos  8 


Figure  1  -  Geometrical  Configuration  of  Cable  System 
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-* 


BT  =0 
WB  =  20.0  lb 

Vbh  -  °-5  ff2 

sv  ■  °-'i  “2 


Buoy  Parameters: 
V?F  =0 

Vs,  •  °-79  f(2 

VlR  "  m-°  lb 


MSV!I 

=  5.19  slugs 

MSVy 

=  6.59  slugs 

RT 

=  RT  =  0 

X 

z 

RT 

=  1.0  ft 

y 

It 

O 

t-h 

rf 

v— <* 

Vol  (ft3) 

CDASH 

AYSM 

AXSH 

AZSM 

1.00 

0.00 

0.00 

0.30 

1.18 

0.746 

0.15 

1.63 

0.746 

0.0 

2.09 

0.746 

-0.15 

2.56 

0.746 

-0.30 

3.01 

0.746 

-1.00 

4.19 

1.00 

Cable  Parameters: 


L  =  100  ft 

C_  = 

0.02 

T  = 

0  lb 

T 

r 

d  =0.12  in 

W 

0.05 

lb/ft 

Cl  = 

2000  lb 

C  =  1  4 

y 

0.02 

si/  f  t 

NCAR  = 

4 

Figure  3  -  Parameters  for  Free-Floating  Buoy-Cable  System 
of  Sample  Problem  1 
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Cable  Parameters: 

L  =  500  ft 
d  =  0.5  in 
CD  -  1.8 

cT  =  0 

Lower  Array: 

L  =  200  ft 
d  =  1.0  in 


W  =  0.2  lb/ft 

y  -  0.0089  sl/ft 

T  =  0  lb 
r 

Cl  =  10000  lb 


W  =  0  Ib/ft 

U  =  0.0109  sl/ft 

T  =  0  lb 
r 

Cl  =  10000  lb 


NCAB  =  5 


Figure  4  -  Parameters  for  Cable  Towing  Lower  Array  of 
Sample  Problem  2 
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Steady-state  calculation  for  continuous  cable 

Dynamic  calculation  for  segmented  cable,  with  and  without  cable  inertia 


TABLE  2  -  VALUES  OF  VELOCITY  x  IN  FT/SEC  CALCULATED  BY  PROGRAM  CAMBOD 
DURING  DEPLOYMENT  OF  FREE-FLOATING  BUOY-CABLE  SYSTEM 
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CAMBOD  ET/RT  =  1.20 
CAMBOD1  ET/RT  =  1.10 

CABMOD°  =  CABMOD  with  no  cable  added  inertia  (AMC 
CABMOD1  =  CABMOD  with  cable  added  inertia  (AMC  = 


TABLE  3  -  (continued) 


Segm 

CAB- 

MOD1 

t  =  1 

CAB¬ 

MOD0 

,500  sec 

CAB- 

3DYN 

CAB¬ 

TRAN 

CAB- 

M0DSS 

1 

39.52 

39.52 

— 

39,10 

39.58 

2 

39.66 

39.66 

— 

38.83 

39.76 

3 

40.16 

40.16 

— 

37.93 

40.37 

4 

47.61 

47.61 

— 

30.05 

44.77 

5 

65.55 

65*56 

— 

71.63 

66.97 

CABMOD1 

ET/RT 

*  0.43 

CABMOD0 

ET/RT 

=*  0.55 

CAB3DYN 

ET/RT 

■  0*63 

CABTRAN 

ET/RT 

»  0.01 

CABMOD1 

=  CABMOD  with  added 

cable  inertia 

CABMOD0 

=  CABMOD  with  no  added  inertia 

ss 

C.ABMOD  =  Steady-state  results  calculated  by  CABMOD 
taking  the  cable  to  be  continuous 
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TABLE  4  -  STEADY-STATE  VALUES  OF  CABLE 

ANGLE  $v  IN  DEG  CALCULATED  BY  PROGRAM 

CABTRAN  FOR  VARIOUS  VALUES  OF  V  AND  NCAB 

S 


NCAB  - 

5 

ET/RT  -  0.01 

Vg  (knots) 

Segm  _ 

1 

3 

5 

10 

1 

39.10 

73.51 

82.76 

87.59 

2 

38.83 

67.53 

73.88 

81.08 

3 

37.93 

79.16 

85.11 

87.99 

4 

30.05 

61.04 

71.91 

80.79 

5 

71.63 

83.83 

86.29 

88.15 

NCAB  - 

10 

ET/RT  •  0.01 

Segm  . 

1 

3 

5 

10 

1 

39.34 

71.64 

78.58 

82.53 

2 

39.32 

71.64 

79.59 

86.71 

3 

39.30 

71.60 

77.72 

81.84 

4 

39.26 

71.83 

80.93 

87.40 

5 

39.19 

70.90 

75.85 

81.16 

6 

39.04 

73.84 

83.34 

88.03 

7 

38.71 

66.86 

73.14 

80.61 

8 

37.61 

80.36 

86.01 

88.48 

9 

27.15 

59.46 

70.92 

80.29 

10 

76.58 

85.51 

87.30 

88.65 
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TABLE  4  (continued) 


NCAB  -  20 

Vg  (knots) 


ET/RT  -  0.04 


Segm 

1 

3 

5 

10 

1 

39.38 

71.64 

78.92 

84.30 

2 

39.38 

71.64 

78.92 

84.65 

3 

39.38 

71.64 

78.92 

84.16 

4 

39.38 

71.64 

78.92 

84.84 

5 

39.38 

71.64 

78.92 

83.91 

6 

39.37 

71.64 

78.92 

85.15 

7 

39.37 

71.64 

78.91 

83.54 

8 

39.36 

71.64 

78.95 

85.60 

9 

39.36 

71.64 

78.84 

83.01 

10 

39.35 

71.64 

79.10 

86.20 

11 

39.33 

71.64 

78.54 

82.35 

12 

39.31 

71.64 

79.66 

86.93 

13 

39.29 

71.60 

77.59 

81.59 

14 

39.24 

71.86 

81.15 

87.68 

15 

39.16 

70.80 

75.53 

80.86 

16 

38.99 

74.10 

83.79 

88.37 

17 

38.61 

66.33 

72.57 

80.25 

18 

37.34 

81.28 

86.69 

88.85 

19 

24.67 

58.25 

"0.17 

79.91 

20 

80.34 

86.77 

88.06 

89.03 
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CABMOD  =  CABMOD  with  no  cable  added  inertia 


TABLE  6  -  VALUES  OF  CABLE  ANGLE  4>v  IN  DEG  DURING 
DEPLOYMENT  OF  CABLE  TOWING  NO  ARRAY  AT  5  KNOTS 
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CABMOD  -  CABMOD  with  no  cable  added  inertia 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMUHANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE,  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


